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1  Background  and  Aim  of  Project 

An  insulation  degradation  of  electrical  cables  of  instruments  and  control  facilities  is  one  of  the 
critical  phenomena  for  aging  management  of  large  scale  systems,  such  as  in  airplane,  power 
plants,  etc  [1] .  From  the  fact  that  dielectric  properties  of  materials  used  in  electrical  cables  can 
be  degraded  at  the  course  of  long-term  service,  use  of  microwaves  includes  potential  applications 
in  nondestructive  test  for  cable  degradation.  The  idea  has  been  already  applied  to  the  in- 
service  inspection  in  the  variety  of  industrial  applications.  Time-Resolved  Microwave  Dielectric 
Absorption  (TRMDA)  enables  to  measure  evanescent  electrical  held  leaking  from  the  probe  into 
the  cable  insulation  by  using  cavity  resonator  [2],  This  electrical  held  changes  with  the  dielectric 
properties  of  material,  which  in  turn  changes  the  resonant  frequency.  Figure  1  depicts  the  overall 
configuration  of  the  current  nondestructive  testing.  The  nondestructive  test  is  performed  by 
evaluating  a  bunch  width  of  rehected  microwave.  This  is  because  the  increasing  the  whole 
of  cavity  in  cable  insulation  can  be  detected  as  variation  of  the  bunch  width.  Although  the 
method  has  been  a  conventional  method  for  characterizing  aging  properties,  there  are  several 
disadvantages  on  the  practical  implementation.  One  is  that  the  detection  and  characterization 
of  variation  of  bunch  width  are  not  feasible  for  on-site  inspection.  Also  continuous  surveillance 
of  the  variation  is  very  crucial  tasks.  Taking  into  remind  these,  our  concern  in  this  paper  is  to 
develop  a  computational  method  for  identifying  dielectric  properties  of  electrical  cables  based 
on  microwave  guide  measurement  system. 

The  quantitative  nondestructive  evaluation  (QNDE)  of  material  aging  parameters  continues 
to  be  a  very  challenging  problem.  In  our  approach,  we  formulate  a  forward  problem  arising 
in  specific  micro  guided  wave  test.  More  specifically,  it  represents  a  mathematical  model  of 
specific  nondestructive  test  using  the  parameter-to-output  mapping  with  the  appropriate  ad¬ 
missible  class  of  material  parameters,  To  this  end,  the  mathematical  model  is  described  by  a 
two-dimensional  Maxwell’s  equations.  A  numerical  scheme  of  the  forward  problem  is  derived 
using  the  finite-difference  time-domain  method  (FDTD)  in  two  spatial  dimensions.  Thus,  the 
inversion  of  material  parameters  is  implemented  with  the  aid  of  the  forward  problem.  The  non¬ 
linear  least  square  identification  (OLSI)  is  a  conventional  inverse  methodology  and  there  have 
been  so  many  efforts  in  variety  of  industrial  applications.  For  OLSI  related  to  electromagnetic 
inversion,  we  refer  [3]  [4]  [5]  [6]  [7].  However, it  is  well  known  that  the  problem  mentioned  above 
has  many  solutions  due  to  the  fact  it  is  ill-posed.  Recently,  interest  has  grown  in  stochastic 
inversion  using  Markov  Chain  Monte  Carlo  (MCMC)  methods  [8]  [9] .  Some  previous  efforts  on 
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minimums  of  the  inverse  solution.  In  this  paper,  a  stochastic  inversion  technique  based  on  the 
Metropolis- Hasting  algorithm  is  applied  to  the  problem  presented  here. 
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Figure  1:  Nondestructive  test  using  TRMDA 


2  Methodical  Procedures 

2.1  Milestone  1:  Formulation  of  nondestructive  evaluation  arising 
in  microwave  analysis 

Guided  microwaves  analysis  is  considered  for  the  detection  and  the  characterization  of  corrosion 
under  insulation.  First,  the  forward  problem  is  formulated  by  a  Maxwell  equation  in  two  or  three 
spatial  dimensions.  Secondly,  the  forward  problem  is  considered  by  constructing  deterministic 
parameter-to-output  mappings  based  on  the  frequency  dependent  media  with  unknown  material 
quantities  as  well  as  with  the  unknown  geometries.  Thirdly  numerical  scheme  is  developed  by 
using  FDTD  method  in  order  to  solve  the  forward  problem. 


2.1.1  Problem  Domain 

(Term  :  From  September  2010  to  December  2010) 


For  the  simplicity  of  discussions,  we  consider  a  two  dimensional  bounded  domain  x  = 
(xi,x2)  G  D  C  R2  as  an  inspection  area.  Let  Dc  be  a  cable  domain  defined  on  fh  The 
cable  domain  consists  of  a  chemical  material  Q,/  and  a  metal  core  fl0.  The  chemical  material 
corresponds  to  a  polymer  of  which  the  dielectric  constant  er  and  conductivity  often  vary 
at  different  frequencies.  Figure  2  depicts  the  two-dimensional  problems  space  considered  here. 
Our  electromagnetic  interrogation  is  then  represented  by  electric  and  magnetic  fields  in  O.  Let 
D3(t,x),  Hi(t,x),  H2(t,x)  and  E3(t,x )  be  the  component  of  electric  flux  density,  magnetic  and 
electric  held  at  time  t  G  [0,  T]  and  at  location  iGfl.  In  order  to  use  Gaussian  units,  we  also 
use  the  normalizing  terms, 


Problem  Domain  Q 


where  eo  and  // o  are  permittivity  and  magnetic  permeability  of  air.  From  Maxwell’s  equations, 
the  electromagnetic  propagation  in  the  transverse  magnetic  (TM)  mode  is  governed  by 
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with  zero  initial  states  and  with  the  absorption  boundary  conditions.  The  complex  relative 
dielectric  constant  e*  of  the  polymer  Qj  can  be  represented  by  a  lossy  dielectric  medium  of  the 
form; 


e 


* 

r 


er  + 


Od 

jue  o ' 


(5) 


On  the  other  hand,  the  value  of  this  parameter  in  Oa  =  0  —  Oc  becomes  e*  —  1.  The  value  of 
that  in  f20  is  also  given  by 
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where  aQ  denotes  the  conductivity  of  the  metal  core.  The  source  current  Js  is  the  test  signal  in 
our  inspection  process  and  is  given  by 

Js(t,  x)  =  8(ns)gs(t)I(0,ts)(t)  (7) 

where  5  denote  the  array  of  actuator  in  such  as  =  0-  In  Eq.  (7),  the  test  signal  is 

truncated  at  a  finite  time  ts  by  the  indicator  function  I.  The  detecting  signal  can  be  made 
through  the  allocation  of  sensor  arrays  h20  in  hi  by 

Y(t,  x;  q)  =  E3(t,  x;  q)  in  HD  (8) 

where  the  parameter  vector  q  to  be  identified  is  denoted  by 

q  =  {er,ad}.  (9) 

Thus  the  direct  problem  for  our  NDT  is  formulated  as 

q^E(f,x;q).  (10) 

2.1.2  Numerical  Scheme 

(Term  :  From  January  2011  to  March  2011) 

The  finite-difference  time-domain  method  (FDTD)  in  two  spatial  dimensions  can  be  im¬ 
plemented  to  construct  the  forward  model  [12].  Putting  Eqs.(l),  (3),  and  (4)  into  the  finite 
differencing  scheme  results  in  the  following  difference  equations: 
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At 

In  order  to  prevent  spurious  reflections  from  the  edge  of  the  problem  space,  the  perfectly 
matched  layer  (PML)  is  adopted  to  two  dimensional  domain  (  See  [13]  for  more  details  ).  The 
principle  of  micro  guided  wave  testing  can  be  represented  by  the  simulation  of  plane  waves. 
Since  the  propagating  plane  wave  must  not  interact  with  the  absorbing  boundary  conditions  like 
PML,  the  problem  space  can  be  divided  up  into  two  regions  such  as  total  fields  and  scattered 
held.  Figure  3  illustrates  how  the  incident  plane  wave  behaves.  To  simulate  a  plane  wave 
interacting  with  a  cable  object,  the  cable  object  according  to  its  electromagnetic  properties 
must  be  specified.  The  simulation  of  a  plane  wave  pulse  hitting  a  dielectric  cylinder  with 
er  =  30  and  a  =  0.3[m/s]  is  demonstrated  in  Fig.  4.  After  the  time  step  T  =  50,  the  pulse  is 
interacting  with  the  cylinder.  At  time  step  T  —  75  some  of  it  passes  through  the  cylinder  while 
the  other  part  of  it  goes  around  that.  It  could  be  also  recognized  that,  at  T  —  100,  the  main 
part  of  the  propagation  is  being  subtracted  out  the  end  of  the  total  held. 
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Figure  3:  Diagram  of  a  plane  wave  hitting  a  cable  object 


2.2  Milestone  2:  Development  of  new  stochastic  inverse  methodol¬ 
ogy 

(Term  :  From  April  2011  to  June  2011) 

With  the  background  knowledge  of  Bayes  formula,  the  full  probability  model  is  specified 
by  the  deterministic  formula  considered  in  Milestone  1  and  the  stochastic  relationship  for  their 
measurement  strategies.  More  specifically,  the  posteriori  density  function  is  provided  with 
likelihood  functional  update  based  on  iterative  procedures  of  measurements  and  with  a  priori 
probability  densities  with  respect  to  the  unknown  quantities.  The  application  of  extreme  value 
statistics  to  corrosion  analysis  is  also  an  important  factor  of  our  approach. 

Our  problem  is  to  identify  dielectric  parameters  related  to  cable  insulation 

q  =  {er,  ad} 

using  the  information  obtained  through  a  scattering  experiment  for  the  forward  problem  stated 
in  the  previous  section.  The  observations  consist  of  sampling  data  associated  with  the  electrical 
field.  We  suppose  that  the  experimental  observations  consist  of  the  values  of  electrical  field  at 
the  measurement  point  xp  =  {x*,Xp}  at  the  sampling  time  interval  At.  Let  Y^  =  {Yd[ri]}^=l 
denote  the  collection  of  the  data  we  seek  to  reconstruct  and  YJV(q)  =  {Y(nAt,xp]  q)}n=i  be 
the  electrical  field  arising  from  the  forward  problem  with  dielectric  parameter  vector  q.  Then 
the  structural  equation  model  is  formulated  as 

Yd[n\  =  Y (nAt,  xp\  q)  +  rjn,  rjn  ~  A/"(0,  r2), 

(n=  1,  2,  •  •  •  N). 
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(a)  T  =  25 


(b)  T  =  50 


(c)  T  =  75 


(d)  T  =  100 


Figure  4:  Simulation  of  a  plane  wave  pulse 


where  r]n  denotes  the  mutually  independent  sequence  of  additive  measurement  noise.  With  the 
background  knowledge  of  Bayes  formula,  the  posteriori  density  function  with  respect  to  the  set 
of  unknown  parameter  vector  q  can  be  simply  written  by 
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where  p(er)  and  p(crd)  denote  a  priori  density  functions  of  er  and  An  estimation  algorithm 
can  be  performed  by  sampling  procedures  for  the  posteriori  distribution  given  by  Eqs.  (12) 
and  (13)  from  which  sample  paths  can  be  drawn  using  Markov  chains.  For  the  practical 
implementation  of  MCMC,  Metropolis- Hasting  algorithm  is  applied  to  the  problem  considered 
here. 


2.3  Milestone  3:  Feasibility  studies  based  on  computational  exper¬ 
iments 

(Term  :  From  July  2011  to  September  2011) 


The  validity  and  feasibility  of  our  approach  developed  in  Milestones  -1  and  -2  will  be  demon¬ 
strated  through  computational  experiments  for  appropriate  specific  examples.  The  Gibbs  sam¬ 
pler  or  the  other  appropriate  MCMC  samplers  (e.g.,  Metropolis-Hasting  algorithm,  hierarchy 
method,  etc)  are  tested  for  the  above  problems.  The  applicability  of  a  real  problem  to  our 
analysis  is  also  considered. 

To  test  the  feasibility  of  the  estimation  approach,  we  produce  synthetic  data  for  the  ob¬ 
servations  YA  =  (T  [n] }  by  adding  random  noise  to  the  results  of  the  simulation  with  a 
known  set  of  parameters  q/rue  =  {e*rue,  The  data  to  be  identified  are  generated  with  the 

forward  problem  in  the  previous  section  using  the  following  parameter  values: 


elrue  =  30.0 
<?Te  =  0.3  [S/m] 

The  carrier  frequency  of  the  plane  wave  is  90 GHz  and  the  duration  of  sampling  time  is  taken 
as  At  =  [s],  respectively.  The  observation  data  consist  of  512  measurements  of  the  electric 
held  taken  at  xp  =  (0.5,  0.5).  Figure  5  depicts  the  model  data  with  true  parameters.  A  priori 
densities  of  er  and  a  are  preassigned  as  normal  distributions; 
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where  e°  and  a0  denote  the  nominal  values  of  dielectric  parameters.  The  artificial  noise  term 
was  provided  using  a  standard  Gaussian  random  generator  iV(0,0.1).  After  10000  iterations 
and  discarding  the  first  1000  iterations,  the  posterior  means  of  q  were  evaluated.  Typical  paths 
based  on  9000  draws  are  presented  in  Fig.  6.  The  resulting  marginal  posteriors  appear  in  Fig.  7. 
Estimation  results  are  summarized  in  Table  1. 


3  Summary  of  Achievement 

Guided  microwaves  analysis  was  considered  for  the  detection  and  the  characterization  of  cable 
insulation.  The  forward  problem  in  two  spatial  problem  domain  was  considered  by  construct¬ 
ing  deterministic  parameter-to-output  mappings  based  on  the  frequency  dependent  media  with 
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Figure  5:  Synthetic  data  with  true  parameter  values  (noise  free) 


Table  1:  Estimation  summary 


Quantities 

6r 

ad[S/m 

True 

30.00 

0.3 

Min 

29.67 

0.2179 

1st  Quater 

31.98 

0.2838 

Median 

32.49 

0.2986 

Mean 

32.51 

0.3000 

3rd  Quarter 

33.04 

0.3154 

Max 

35.49 

0.4157 

material  quantities  of  cable  insulation.  Numerical  scheme  is  developed  by  using  FDTD  method 
in  order  to  solve  the  forward  problem.  Secondly,  stochastic  inverse  methodology  was  considered 
with  the  background  knowledge  of  Bayes  formula.  The  full  probability  model  was  specified  by 
the  deterministic  formula  treated  in  the  forward  problem  and  by  the  stochastic  relationship 
for  their  measurement  strategies.  The  posteriori  density  function  was  provided  with  likelihood 
functional  update  based  on  iterative  procedures  of  measurements  and  with  a  priori  probability 
densities  with  respect  to  the  unknown  quantities.  The  validity  and  feasibility  of  our  proposed 
method  were  demonstrated  through  computational  experiments  for  appropriate  specific  ex¬ 
amples.  The  MCMC  samplers  using  Metropolis-Hasting  algorithm  was  tested  for  our  inverse 
problem. 

4  Refereed  Publications 

Fumio  Kojima  and  J.S.  Knopp,  “Quantitative  Evaluation  for  Two-dimensional  Electromagnetic 
Propagation  in  a  Dielectric  Medium  using  Markov  Chain  Monte  Carlo  Method  ”  ,  Proceedings 
of  the  43th  ISCIE  International  Symposium  on  Stochastic  System  Theory  (SSS  ’  11)  ,  Otsu, 
Japan,  October  27-28,  2011  (Submitted). 
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Figure  6:  Trajectory  of  the  chains 


5  Collaborations  related  in  AOARD  Work 


[1]  Professor  H.  T.  Banks,  Director,  Center  for  Research  in  Scientific  Computation,  North 
Carolina  State  University,  Raleigh  NC  USA 

Research  collaboration  has  been  performed  during  April  17-21,  2011.  Project  on  mate¬ 
rial  damage  parameter  estimation  for  electromagnetic  acoustic  interaction  has  been  started. 
We  will  accomplish  the  convergence  properties  of  inverse  solution  of  nondestructive  test  using 
electromagnetic  acoustic  transducer. 

[2]  Mr.  Jeremy  S.  Knopp,  Researcher,  AFRL/RXLP,  Nondestructive  Evaluation  Branch,  Ma¬ 
terials  and  Manufacturing  Directorate,  Air  Force  Research  Laboratory,  Dayton  OH  LISA 

Research  collaboration  has  been  performed  during  April  22-25,  2011  in  Dayton.  Also  he 
visited  Kobe  University  during  January  26-  28,  2011.  We  discussed  new  inverse  methodologies 
for  structural  health  monitoring  which  includes  how  statistical  approach  by  merging  modeling 
of  material  failure  process  and  how  to  design  inverse  methodologies. 
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